Integrable Mott insulators driven by finite electric field 
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We develop a method for extracting the steady nonequilibrium current from studies of driven 
isolated systems, applying it to the model of one-dimensional Mott insulator at high temperatures. 
While in the nonintegrable model the nonequilibrium conditions can be accounted by internal heat- 
ing, the integrability leads to a strongly nonlinear dc response with vanishingly small dc conductivity 
in the linear-response regime. The finding is consistent with equilibrium results for dc limit of the 
optical conductivity determined in the presence of a weak and decreasing perturbation. 

PACS numbers: 72.10.-d, 71.27. +a,72.10.Bg 



Introduction. — Understanding the differences between 
integrable and nonintegrable systems has been recognized 
as a mostly theoretical problem. However, recent ad- 
vances in the physics of ultracold atoms enable direct 
experimental verification of the underlying theoretical 
concepts. The presence of conserved quantities (a hall- 
mark of integrability) restricts thermalization processes 
[H , hence integrable systems may relax to the generalized 
Gibbs state [2| instead of the thermal state. Qualitative 
differences between both system types show up also in 
their linear-response (LR) to an external field F. Al- 
though integrability is usually broken by F, this breaking 
is not visible in the LR theory, since LR susceptibility is 
calculated at F = 0. The role of integrability at finite 
driving remains important but largely unexplored field. 
Here we study insulators driven by electric field F and 
investigate how the integrability-related properties of the 
LR regime decay/change with increasing driving. 

In the integrable metallic systems there is a dissipa- 
tionless component of the LR optical conductivity er(cj) 
at arbitrary temperature T 0]. The difference between 
integrable and nonintegrable metals persists also beyond 
the LR regime but only for restricted time of driving 0] ■ 
The response of integrable systems is dominated by the 
Bloch oscillations which, however, are damped due to 
field-induced breaking of integrability. 

In case of Mott insulators at T > 0, LR has proven 
to represent by itself a theoretical challenge not fully un- 
derstood so far, with a transport fundamentally different 
from a usual picture of a band insulator at T > 0. It is 
plausible that generic nonintegrable models do not show 
dissipationlcss transport at T > as characterized by the 
charge stiffness D(T) > 0. This seems to be quite settled 
also for some prototype one-dimensional integrable mod- 
els within the insulating regime [1, 0], as e.g. the t—V (or 
anisotropic Heisenbcrg) model. However, in the latter 
case also the regular dynamical conductivity a{uj — > 0) 
vanishes in finite systems at T > Q (stimulating the 
scenario of an 'ideal insulator' Q), but a(oj) shows at the 
same time anomalously large finite-size contribution at 
u) oc 1/L which at the first sight might be consistent with 



thermodynamic value cr(0) > 0. It should be noted that 
tr(0) > may appear as a result of steady-current calcu- 
lations for an insulator with attached leads as well as 
of the time-evolution analysis Q . In view of these open 
questions , it is quite unclear which might be the ef- 
fect of finite but weak driving F > as well as of small 
nonintegrable perturbation. 

In this Letter we show that for nonintegrable insula- 
tor one can indeed reconcile the response to finite-field 
driving with the LR result if we take into account the 
energy increase (the Joule heating) and the existence 
of (time-dependent) quasi-cquilibrium after a short-time 
transient. More dramatic are results for driven integrable 
insulators. Since F ^ breaks integrability large F- 
response is quite similar to the nonintegrable case. On 
the other hand, our results for weak field F — > indicate 
the vanishing of the steady current response cr(0) = I/F 
reviving the claim of an 'ideal insulator' with ct(0) — > 
at T > 0. Such scenario seems to follow also from the al- 
ternative approach calculating LR er(0) with a weak and 
decreasing integrability-breaking perturbation. 

Model. — A one-dimensional system of interacting 
spinlcss fcrmions with periodic boundary conditions is 
threaded by time-dependent magnetic flux cj){t) 

H = -^{e^c^+h.c.} 
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+V H "A+l + W (1) 
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where n,- 



c\cj 



th is the nearest-neighbor (nn) hopping 
integral, V and W are the nn and next nn repulsive po- 
tentials, respectively In the absence of external driving 
the system is integrable for W — 0, whereas nonvanishing 
W breaks the integrability. As we investigate the physics 
of Mott insulators, we consider a half-filled system with 
L sites and L/2 fermions. Initially <f>{t = 0) = and 
the system is assumed to be in an (approximate) micro- 
canonical state. Then, for t > the flux <fi(t) = —Ft 
gives rise to a constant electric field F. The field will be 
expressed in units of [th/eoa], where eo is the unit charge 
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FIG. 1: (Color online) Nonintegrable case V = 3, W = 1, 
L — 26 with a pure initial state (N a = 1). (a): current I(t) 
(main) and I(t)/F (inset) for initial /3 e ff ~ 0.42; (b) and (c): 
I/F vs. fteff for various initial energies Eo (panel b) and 
various F (panel c); (d): The same as in c) but for various L 
and compared with the results of equilibrium LR theory. 



and a is the lattice distance. We set th = ft = eo = a = 1. 
The charge current I(t) can be calculated cither directly 
from the expectation value of the current operator: 



J 



L 



£{e^>ct +1 c,-h.a} 
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or from the increase of the total energy Q 



E(t) = -(H(t))=LF(t) 1(f). 



(2) 



(3) 



Calculating the dc current / for finite electric field is 
a challenging problem. Literally, the dc response (steady 
nonequilibrium state) is expected only for open quantum 
systems when the energy gained by carries moving under 
F [see Eq. ([3])] is transferred to the heat reservoirs. In 
driven isolated systems however, the energy steadily in- 
creases causing a steady decay of I(t). This holds true 
independently of intcgrability or the system size and re- 
mains valid also in the thermodynamic limit. Wc first 
analyze nonintegrable systems and show how I can be 
reliably determined from the decaying I(t). 

Nonintegrable case. — As the technical aspects of 
our computational approach are described in Ref. 
here we recall only the basic steps. We apply micro- 
canonical Lanczos method [l0 ] and numerically gener- 
ate N s initial states ^(O)) with assumed energy E n = 
(H(t = 0)) but as well with a small energy uncertainty 
SE = {[H(0) - Eq] 2 ) 1 /' 2 ~ 10" 2 . Then, the time evolu- 
tion of l'I'z(t)) is calculated by step-vise change of the 
flux <j)(t) in small time increments St 1 applying the 
Lanczos time propagation method 12( to each interval 



St. The larger the system is, the fewer initial states are 
needed. It will be shown that for large L one can take 
N s = 1 and omit the averaging over |y;(0)). Similar ob- 
servation has been reported in Ref. |13j but for relaxation 
instead of evolution upon driving. 

Inset in Fig. [1^ shows the ratio I(t)/F for a nonin- 
tegrable system. After a relatively short time of driving 
(t ~ 2) results for different fields split marking the end 
of the equilibrium LR regime. This nonlinearity becomes 
even more important for longer driving, when the largest 
currents occur for the weakest fields (see main panel). In 
the following we show that such behavior occurs simply 
due to heating. For this sake it is convenient to introduce 
a renormalized energy e(t) = [Eoc — E(t)]/L, where 



L(V + W) 



L/2- 1 



(4) 



denotes the energy for T — > oo. In equilibrium and at 
high temperatures e is proportional to the inverse tem- 
perature /3 



v 2 + w 2 

16 



(5) 



We use an effective inverse temperature f3 e f / defined by 
Eq. (O also for driven systems, when this quantity 
should be understood only as a measure of the instan- 
taneous energy. Upon driving f3 e ff(t) monotonically de- 
creases provided I(t) does not change sign. 

Fig. QJ shows the same data as panel (a) but the ra- 
tio I/F is now plotted as a function of /? e //(i) instead 
of t. After the system passes the transient state with a 
characteristic peak in I(t), results for different F nearly 
overlap. Hence, the proportionality between / and F is 
restored far beyond the LR regime provided that the ef- 
fect of the Joule heating is appropriately subtracted. It 
also explains the origin of nonlinearity visible in Fig. QJi: 
fixing the time of driving while allowing for different F 
one actually compares systems with exceedingly different 
energies. Such inconsistency is not a finite-size effect but 
is generic for any isolated system with a non-vanishing 
carrier density. As shown in Fig. [TJ), I((3 e ff)/F is al- 
most independent of the initial energy E or the choice 
of \^>i (0)). This ratio is uniquely determined by the in- 
stantaneous energy (f3 e f f) up to fluctuations which van- 
ish in the thermodynamic limit (see panel d). Fig. [TJi 
also demonstrates that 7[/3 e //(i)] is close to the the dc 
current obtained from the equilibrium LR as described 
in more detail later. 

Result presented in Fig. Q] consistently support the fol- 
lowing picture: the isolated system evolves upon finite F 
towards the infinite-temperature state but the (quasi) dc 
current is uniquely determined by the instantaneous en- 
ergy, suggesting a quasi-equilibrium evolution. Attach- 
ing the system to a heat reservoir should set /3 e / / leading 
to I[P e ff(t)] = const. However, in such case / could de- 
pend on the details of the coupling to reservoir. This 
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FIG. 2: (Color online) (a) and (b): I/F vs. /3 e // for integrable 
system V = 3, W = and L = 28. (c) and (d) /3 e// vs. tF 2 
for integrable (c) and nonintegrable (d) systems. In (c) curves 
are shifted horizontally so they overlap for large times. N s = 1 
apart from F = 0.15 and F = 0.2 in (a) where N s — 7. 
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FIG. 4: (Color online) Equilibrium LR: a(cu)/P in the limiting 
case ft — > for various IV" (panels a and c) and various system 
sizes L (panel b). Main panel in d) shows <r(0)//3 vs. 1/L 
(points) while inset shows W-dependence of cr(0)//3 for V = 3 
extrapolated to L — > oo . 
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FIG. 3: (Color online) I/(/3 eff F) estimated from E(t) [Eq. 
©] for fieff — > (points with guidelines). For nonintegrable 
case results for /3 e // =0.1 (diamonds) together with LR for 
13 = 0.1 (dotted line) are shown as well. 



observation allows us to estimate intrinsic I for finite F 
even though calculations are carried out on an isolated 
system. As expected for nonintegrable insulators, the 
Joule heating represents the central source of nonlinear- 
ity at high temperatures. 

Integrable case. — Before applying the same approach 
to the integrable insulators one should address the ques- 
tion whether / can be again uniquely determined (up 
to finite-size fluctuations) by the instantaneous energy. 
Since finite F breaks integrability such scenario is possi- 
ble and as follows from Fig. [2Jd it actually takes place. 
Although, for a pure initial state (N s = 1) the finite 
size fluctuations are much larger than in the noninte- 



grable case (Fig. [TJ;), I((3 e ff) for various initial states 
do overlap. Therefore, one can apply the same method 
of reasoning as for nonintegrable insulators and take into 
account the heating (the energy increase). The results 
are shown in Fig. [2ji. Here, one finds important qualita- 
tive differences between both system types. While in the 
nonintegrable case the ratio I/F slightly decreases with 
F (see Fig. QJ), now an opposite dependence is visible. 
Namely, I((3 e ff)/F (strongly) increases with F indicat- 
ing on a nonlinear mechanism beyond that of the Joule 
heating. If the dependence between F and / were linear 
[I(E, F)/F = -/(E)], one would obtain from Eq. © that 
dE/dt oc F 2 "f(E) and all curves in Fig. would overlap. 
Unlike results for nonintegrable insulators shown in Fig. 
[2}i, for integrable case the slope of E(t) increases faster 
than F 2 . Although the accessible values of L do not al- 
low for a detailed size scaling, Fig. [2t strongly suggests 
that the nonlincarity is not a finite-size effect. 

Fig. [3] shows our central result: F dependence of the 
ratio I/F for /3 e ff <C 1. In the case of nonintegrable in- 
sulators, I/F approaches LR a(Q) when F -> (dotted 
horizontal line) , while weak but systematic decreasing of 
this ratio occurs for larger F. This is consistent with a 
nonmonotonic I—F characteristic which has been demon- 
strated for various tight-binding models 



14 . However, 



in integrable insulator the ratio I/F depends on F in a 
qualitatively different way and LR regime is not visible 
down to the lowest accessible fields. 

Since the difference due to integrability is most pro- 
nounced for low electric fields, the underlying mechanism 
should be visible also in the LR. Fig. @]shows equilibrium 
optical conductivity obtained within the microcanonical 
Lanczos method introduced in Ref. 0, EH ■ As already 
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noticed for the integrable insulator (V > 2) Q tr(o; > 0) 
is dominated by a finite-size dependent low frequency 
peak at loq ~ 16/L, its position only weakly varying 
with V. But the existence and the value of <r(0) (at 
any L) critically depends on W (see Figs. andQJi). 
Since W ^ (as well as F > 0) breaks integrability 
one can expect finite tx(0) in the thermodynamic limit. 
However the L — > oo scaling remains delicate as shown 
in Fig. HJd for W = 0.5. For u > lo fs ~ 1.5 all the 
data form a single L-independent curve which could sug- 
gest with an oversimplified extrapolation to L — > oo that 
tr(0)//3 > 0.5. However, the high-T sum rule (at half- 
filling), J Q doja(uj) = which is moreover V,W in- 
dependent, requires L-independent Jq FS dwaiuj), hence 
prevents an extrapolation to e.g. cr(0)//3 > 0.4. A direct 
scaling of a(0)//3 suggests cj(0)//3 ~ 0.15 for W = 0.5 as 
shown in Fig. [4jl. Analogous scaling for other W > 
leads to values of er(0) decreasing and (linearly) vanishing 
with W (see inset in Fig. |4ji). In any case, the result- 
ing tr(0) is much smaller than straightforward scaling of 
anomalous fluctuations for W = which would lead, e.g. 
for V = 3, to (7(0) > 0.3 0. Clearly, the surprising result 
ct(0)(w/->o) — is consistent with the field dependence 
I/F for F — s- shown in Fig. [3J In like man- 
ner as W , finite F as well breaks integrability, leading 
effectively to an enhancement of the response function. 
Consequently / increases with F faster than linearly. 

Discussion. — While our numerical results were obtained 
only for very high T 3> it is plausible that their impli- 
cations might be extended to all T > 0, especially since 
at low temperatures strongly nonlinear I—F characteris- 
tic is well established [l5j]. The notion of 'ideal insulator' 
means vanishing of the LR cr(0) at arbitrary temperature 
that implies at the same time a nonlinear dc response of 
integrable insulators down to the lowest electric fields. 
However, for realistic (open) systems it means that the 
LR regime is determined by the mechanisms which break 
integrability (see discussion in Ref. [1]), hence this regime 
may terminate at low fields. Here, we present two inde- 
pendent results which support the hypothesis of 'ideal 
insulator': nonlinear real-time current for finite driving 
and LR <t(uj) obtained for small but finite integrability- 
breaking parameter W . The Anderson insulator (non- 
interacting case when all states are localized) has been 
the only well known example of an 'ideal insulator'. An 
emerging question concerns the differences between the 
present study and those which suggest <r(0) > A 
possible explanation is that in our setup driving doesn't 
destroy the homogeneous distribution of charge carriers, 
what may happen, e.g., in an open system driven by a 
difference in chemical potentials. We have also studied 
an isolated system exactly at half-filling (microcanonical 
ensemble). Consequently, the Mott insulating phase is 
not destroyed (even locally) by a departure from half- 



filling and integrability is not broken by external leads. 

Our results confirm that the Mott insulators, in partic- 
ular the integrable ones, behave under driving fields qual- 
itatively different from usual band insulators. While in 
the latter case thermally (as well as photo etc.) induced 
carriers, i.e. electrons and holes, are quite independent 
but oppositely charged quasipartilces, in Mott insulators 
the corresponding doublons and holons are highly corre- 
lated, and moreover cannot pass each other in the inte- 
grable ID system (evident in the limit V » 2). Hence 
no dc current can be induced at weak F at any T. On 
the one hand, anomalous low- frequency peak at uj cx 1/ L 
clearly requires caution in the interpretation. On the 
other hand, the concept of the mean free path A would 
be ill defined if the 1 /L dependence in dc response per- 
sists also for L 3> A. 
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